Matter-wave gap solitons in atomic bandgap structures 
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We demonstrate that a Bose-Einstein condensate in an optical lattice forms a reconfigurablc 
matter-wave structure with a band-gap spectrum, which resembles a nonlinear photonic crystal for 
light waves. We study in detail the case of a two-dimensional square optical lattice and show that 
this atomic bandgap structure allows nonlinear localization of atomic Bloch waves in the form of 
two-dimensional matter-wave gap solitons. 
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Recent experiments have demonstrated that 

a Bose-Einstein condensate (BEC) loaded into an optical 
lattice is a perfect test ground for a range of fascinating 
physical effects. In particular, this is because the effective 
potential of an optical lattice can be easily manipulated 
by changing the geometry, polarization, phase, or inten- 
sity of the laser beams. Due to inherent nonlinearity of 
the coherent matter waves introduced by inter-atomic in- 
teractions, BEC in a lattice potential can form a periodic 
nonlinear system, which is expected to display rich and 
complex dynamics. 

On the other hand, propagation of light in dielectric 
structures with a periodic variation of the refractive in- 
dex is receiving growing attention. Photonic bandgap 
(PBG) materials ja] - artificial periodic structures with 
a high index contrast - can be used to effectively control 
the flow of light. The study of photonic crystals made of 
a Kerr nonlinear material, the so-called nonlinear pho- 
tonic crystals, has revealed that such structures exhibit 
a wealth of nonlinear optical phenomena and, in partic- 
ular, they can support self-trapped nonlinear localized 
modes of the electromagnetic field in the form of opti- 
cal gap solitons 0, Q • Dynamically reconfigurable PBG 
structures - optically-induced refractive index gratings in 
nonlinear materials - are now offering new ways to control 
light propagation and localization 

As demonstrated in this letter, BEC in an optical lat- 
tice can be regarded as a fully reconfigurable analog of a 
nonlinear photonic crystal for matter waves - an " atomic 
band gap" (ABG) structure. A deep analogy between co- 
herent light and matter waves suggests that the concepts 
employed in the study of nonlinear PBG structures can 
be borrowed for analyzing BEC in optical lattices. Many 
properties of ABG structures can potentially be exploited 
for high-precision control and manipulation of coherent 
matter waves in a similar way to how PBG structures are 
used to manipulate light. The effect of nonlinear local- 
ization in bandgaps is one of these properties. 

Localized nonlinear excitations of coherent matter 
waves - bright atomic solitons - could be very useful 
for applications such as atomic interferometry due to 



their robust nature. However, their generation has so 
far been experimentally achieved only in BECs with at- 
tractive atomic interactions (analogous to self-focusing 
optical media) which are unstable against collapse above 
a small (~ 10 3 ) critical number of particles. 

In theory, a shallow ID optical lattice can support 
bright matter-wave solitons even in repulsive BEC [t| 
with large atom numbers. These solitons, described in a 
framework a coupled- mode theory jjj, are localized on a 
large number of lattice wells, and are predicted to exist 
only in atomic bandgaps. In the opposite case of a deep 
lattice, the condensate can be described by the superpo- 
sition of ground state modes in the individual wells. The 
mean-field treatment of the condensate in this regime 
leads to a discrete equation which admits solutions in 
the form of stationary modes localized on a few lattice 
sites - discrete solitons |lfjllll| . in complete parallel with 
spatial solitons in periodic optical structures 01 an d lo- 
calized modes of atomic lattices 

The theory of the nonlinear localized matter waves in 
optical lattices is mostly limited to ID case [2jj. How- 
ever, as can be deducted from the analogous studies of 
2D PBG structures [y], the nonlinear localization of in 
BEC in higher-dimensional lattices is qualitatively dif- 
ferent because both the symmetry and dimensionality of 
the lattice start to play an important role in the forma- 
tion and properties of the bandgap structure and cor- 
responding nonlinear modes. In this Letter, we analyze 
the as yet unexplored problem of the existence and sta- 
bility of 2D matter-wave gap solitons of BEC with repul- 
sive interatomic interactions loaded into optical lattices, 
which are analogous to localized states of light waves in 
nonlinear photonic crystals. We show that their accu- 
rate description is only possible within a full mean-field 
model of BEC in a periodic lattice potential, beyond the 
tight-binding approximation or coupled-mode theory of- 
ten employed for the study of BEC in ID optical lattices. 

The dynamics of a three-dimensional BEC cloud 
loaded into a two-dimensional optical lattice can be de- 
scribed by the Gross-Pitaevskii (GP) equation for the 
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macroscopic wavefunction, *f?(r, t) of the condensate, 
.5* 



dt 



1. 



V(r)+ l3D \y\< 



(1) 



This equation has been made dimensionless by using the 
natural length, time, and energy scales of the optical 
lattice. The characteristic length scale in this model, 
aL = k^ 1 = d/n, defines the size of a single well of 
the lattice potential through the lattice constant that de- 
pends on the laser wavelength as d = A/ 2. The depth of 
the lattice potential, measured in the lattice recoil en- 
ergy El = H 2 /(2ma L ), is proportional to the laser inten- 
sity and can be varied from to ~ 20El in experiment. 
The characteristic frequency scale is lol = El/K. The 
interatomic interactions are characterized by the coeffi- 
cient 73D = 87r(a s /az,), where a s is the s-wave scattering 
length, positive for repulsive interactions. 
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FIG. 1: Dispersion diagram for atomic Bloch waves in a 2D 
lattice in the reduced zone representation: shaded - first three 
energy bands, numbered - band gaps (Vo = 5.0). Dashed - the 
line fi = Vo- Right top: contour plot of the lattice potential 
in Cartesian space, black shading corresponds to potential 
minima. Right bottom - the first Brillouin zone of the 2D 
lattice in the reciprocal lattice space, marked are the high- 
symmetry points of the irreducible zone. 

In the simplest case of a square optical lattice, the 
combined potential of the lattice and magnetic trap, 
V(r) = V~L{r) + Vr{r), can be written as V(r) = 



Vo(si 



sm 2 x 



sin y) +oji 2 r 2 /2, where Vq is the amplitude 



of the optical lattice, and w, are the (normalized) trap- 
ping frequencies in the corresponding directions. The 
current experiments with 2D optical lattices are per- 
formed with a strongly anisotropic BEC cloud, with the 
ratio of the trapping strengths fl = u x /oj ytZ ~ 1CP 1 
2J. To simplify our analysis, we assume that the rela- 
tively weak magnetic confinement characterized by trap 
frequencies oji has small effect on the stationary states of 
the condensate in the lattice. Under this assumption, the 
trap component of the confining potential in the lattice 
plane can be neglected, and the model can be reduced 
to a two-dimensional GP equation by assuming separa- 
ble solutions of the form, ip(r;t) = &(z)ip(x,y;t) where 




FIG. 2: Band-gap structure of the atomic Bloch-wave spec- 
trum as a function of the lattice amplitude Vo. Shaded - 
Bloch-wave bands, numbered - gaps. Dashed - the line fi = Vo. 



the function &(z), describing the condensate in the di- 
rection perpendicular to the lattice, is a solution of the 
ID quantum harmonic oscillator problem, with the nor- 
malization condition J $> 2 (z)dz = 1. Then the standard 
dimensionality reduction procedure , leads to the fol- 
lowing equation for the 2D condensate wavefunction in 
the 2D lattice potential: 
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-yl + v L (r) + \M 2 



(2) 



where Vj_ = d 2 /dx 2 + d 2 /dy 2 , Vh(r) = Vh(x,y) is the 
periodic potential of the optical lattice, and the wave- 
function is rescaled as tp — > ip ^/72D , with 72D = 73D / \/2- 
Spectrum of atomic Bloch waves in the optical lattice 
can be found by exploiting the analogies with the theory 
of single-electron states in crystalline solids. Stationary 
states of the condensate in an infinite periodic poten- 
tial of a 2D optical lattice are described by solutions of 
Eq. @ of the form: ip(r,t) = (j>(r) exp(— ifit), where 
\x is the chemical potential, corresponding to the energy 
level of the stationary state in the lattice potential. The 
case of noninteracting condensate formally corresponds 
to Eq. (J2J being linear in ip. According to the Bloch the- 
orem, the stationary wavefunction can then be sought in 
the form </>(r) = Ufc(r) exp(ifcr), where the wavevector 
k belongs to a Brillouin zone of the square lattice, and 
Wfc(r) = Uk{r + d) is a periodic (Bloch) function with 
the periodicity of the lattice. For the values of k within 
an n-th Brillouin zone, the dispersion relation for the 
2D Bloch waves, n n (k), is found by solving the following 
linear eigenvalue problem 



(-iV ± + k) 2 + V L (r) 



(3) 



The eigenvalue problem Q is simpler than a general 2D 
problem due to the separability of the lattice potential. 
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The typical dispersion relation for the energy of the Bloch 
states in the lowest Brillouin zones is shown in Fig. ^ 
for a moderate value of the lattice depth Vq = 5. The 
dispersion diagram is presented in the reciprocal lattice 
space, and the dispersion relations are calculated along 
the characteristic high-symmetry directions of the irre- 
ducible Brillouin zone (see Fig. ^ bottom right). It can 
be seen that the absolute bands and gaps are determined 
by the spectra of the Bloch waves in the middle (point 
r) and on the edge (point M) of the Brillouin zone. 

To understand how the variation in the lattice param- 
eters affects the band structure, we note that, owing to 
the scaling properties of the model, any increase in the 
well spacing, at, translates into decreasing the well depth 
Vo ■ Thus the global behaviour of the band structure can 
be understood by examining the Bloch wave spectrum 
as a function of the lattice depth. Figure [21 presents the 
chemical potentials corresponding to the atomic Bloch 
waves at the edges of the absolute bands plotted as func- 
tions of Vo- The dashed line ji = Vo (dashed line in Fig. 
^) separates the quasi- unbound (fj, > Vo) and strongly 
bound (n < Vo) condensate states in the lattice (see also 
Ref. 15]). By varying the intcrwcll separation or the 
amplitude of the lattice, different regimes of the atomic 
bandgap structure can be accessed. The tight-binding 
regime of the condensate dynamics corresponds to the 
domain Vq V * = /i, where the bands "collapse" to 
discrete levels of bound states in a single isolated well 
(see Fig. EJ). The narrow-gap regime of a shallow lattice, 
where the coupled-mode theory applies, is found for the 
opposite case, Vq <C V *. 

The study of the ID atomic bandgap structures [TiUlfij ] 
indicates that the correct structure and dynamical prop- 
erties of the nonlinear localized modes can only be re- 
vealed through the analysis of the full continuous mean- 
field GP model that bridges the gap between the coupled- 
mode theory and discrete tight-binding approach. One 
promising development in approximate analysis would 
be generalization of the Wannier function formalism |l7^ 
to the 2D case. However, the complexity involved in 
the generation of optimally localized higher-dimensional 
Wannier functions may not provide significant ad- 
vantage compared with the direct solution of the model 
equation J2J. 

Here, we find spatially localized stationary solutions of 
Eq. (0) numerically. Our numerical procedure involves 
minimization of the functional N = J pfdr, where 
f(<t>) = [Vi + n - V L (r) - \</>\ 2 ]<t>, by following a de- 
scent technique with Sobolev preconditioning 19]. The 
minimization procedure yields a stationary state when 
N(4>) — 0. Previous applications of this method to the 
analysis of optical solitons has shown that it is suitable 
for tackling both the fundamental and higher-order non- 
linear localized modes. Our numerical method is inde- 
pendent of the lattice symmetry; it offers a versatile tech- 
nique for finding and analyzing nonlinear localized modes 




FIG. 3: Peak amplitude, <j>a = 0(0,0), of the fundamental 
gap soliton in the first gap (Vo = 5.0); shaded - the first band, 
dashed - line fi = Vo- Below: Spatial structure of the gap soli- 
tons (a) near the high-energy edge of the first band, fi = 3.5, 
and (b) at the top of the lattice potential, /i = 5.0. Contour 
plots on the 3D plots (left) show the structure of the lattice 
potential, darker contours correspond to potential minima. 
Potential contours on the right are an eyeguide. Insets: the 
cross-sections of the localized modes and the potential along 
the line (x,0). 



in a range of different bandgap structures, and for mak- 
ing accurate predictions about their density, degree of 
localization, and spatial structure in different areas of 
the parameter space. 

We have identified different families of matter-wave 
gap solitons of the repulsive BEC in the 2D atomic 
bandgap structure. These solitons can exist in all gaps, 
excluding the semi-infinite gap of the spectrum (marked 
by 1 in Figs ^ 12) below the first band. We note in 
passing that only conventional nonlinear localization of 
BEC with attractive interactions is possible in the semi- 
infinite gap. The change in the peak amplitude of the 
the lowest-order "fundamental" gap solitons across the 
first gap, along with their spatial localization properties, 
are demonstrated in Fig. |3J for Vo — 5. Gap solitons of 
repulsive BEC branch off the lower-energy (i.e. smaller 
values of fx) gap edge (see Figs. |3 top). Near this edge, 
i.e. at low particle numbers and amplitudes, their cen- 
tral peak is localized around a single well [Figs. El (a)]. 
Further within the gap, the localized mode of the BEC 
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contains larger numbers of atoms, and has a strongly 
localized central peak. Near the linear tunnelling thresh- 
old, fj, ~ Vq, the nonlinear localized states develop ex- 
tended tails in the orthogonal directions (x,y), along 
which the tunnelling is assisted by the lower inter-well 
barrier heights (compared to 2Vb in diagonal directions). 
Although the stability analysis of the states localized on 
the extended background (for \i > Vq) is difficult due to 
the large computational domain, we have confirmed that 
strongly localized 2D modes are dynamically stable. In 
the lattice geometry under consideration, the line fj, = Vq 
lies within the first gap (marked 2 in Figs. HIUj an( i only 
localized modes within this gap have been analyzed. 

The higher-amplitude gap solitons have a nontrivial 
structure of the tails [see Fig. (b)], i.e. the zeros of 
the matter wave are centered within individual lattice 
wells (see Figs. [3J insets), rather than on the poten- 
tial maxima. This structure is determined by the spatial 
structure of the Bloch wave at the lower-energy edge of 
the second band (point M in Fig. and cannot be 
described within the framework of the discrete model. 
We have also found families of higher-order gap modes 




FIG. 4: Spatial structure of the higher-order gap solitons 
formed by two (left) out-of-phase, and (right) in-phase fun- 
damental modes (Vo = 5.0, /i = 4.0). Filled contours show 
the lattice potential; dark shading corresponds to potential 
minima. 

that can be identified as bound states of the fundamen- 
tal gap solitons. These states can be centered on the 
lattice potential minima or maxima, similarly to higher- 
order "odd" and "even" states in ID lattices ^ij, and 
can exhibit symmetry-breaking instabilities. Two exam- 
ples of the lowest-order in- and out-of-phase odd modes, 
centered on a lattice maximum, are shown in Fig. 01 

The crucial issue of the potential observation of 2D 
gap solitons is the stability of the 2D localized state 
in a 3D BEC cloud, i.e. dynamical stability of the 
3D state ^(x,y,z;t) with the initial condition given by 
<&(z)(f)(x, y), where <fi(x, y) is the two-dimensional station- 
ary gap soliton. The analysis of stability of 2D solitons 
in a 3D BEC cloud, wich will provide clues for possible 



experimental observation of multi-dimensional gap soli- 
tons, is the subject of our separate study. 

In conclusion, we have demonstrated that the interac- 
tion of BEC with the lattice potential is analogous to the 
light scattering by a nonlinear photonic bandgap struc- 
ture. We have studied the properties of two-dimensional 
atomic bandgap structures and demonstrated the exis- 
tence of gap solitons, the spatially localized states of BEC 
existing in the gaps of the matter- wave spectrum. We be- 
lieve the analogy between the physics of BEC in optical 
lattices and photonic crystals can be useful for reveal- 
ing many novel features of the matter-wave dynamics in 
reconfigurable atomic bandgap structures. 
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